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' — ' ABSTRACT 

^~* HESS J1858+020 is a TeV gamma-ray source that was reported not to have 

^^ any clear cataloged counterpart at any wavelength. However, it has been re- 

^^ cently proposed that this source is indirectly associated with the radio source, 

re-identified as a supernova remnant (SNR), G35.6-0.4. The latter is found to 
!>■ 
f^ be middle-aged (~ 30 kyr) and to have nearby molecular clouds (MCs). HESS 

,_l J1858+020 was proposed to be the result of the interaction of protons accel- 

^ erated in the SNR shell with target ions residing in the clouds. The Fermi 

jC Large Area Telescope (LAT) First Source Catalog does not list any source 

Cd coincident with the position of HESS J1858-I-020, but some lie close. Here, we 

analyse more than 2 years of data obtained with the Fermi-LAT for the region 
of interest, and consider whether it is indeed possible that the closest LAT 
source, IFGL J1857.1+0212c, is related to HESS J1858-I-020. We conclude it is 
not, and we impose upper limits on the GeV emission originating from HESS 
J1858+020. Using a simplified 3D model for the cosmic-ray propagation out 
from the shell of the SNR, we consider whether the interaction between SNR 
G35.6-0.4 and the MCs nearby could give rise to the TeV emission of HESS 
J1858+020 without producing a GeV counterpart. If so, the pair of SNR/TeV 
source with no GeV detection would be reminiscent of other similarly-aged 
SNRs, such as some of the TeV hotspots near W28, for which cosmic-ray diffu- 
sion may be used to explain their multi-frequency phenomenology. However, 
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for HESS J1858+020, we found that although the phase space in principle 
allows for such GeV-TeV non-correlation to appear, usual and/or obscrva- 
tionally constrained values of the parameters (e.g., diffusion coefhcients and 
cloud-SNR likely distances) would disfavor it. 
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1 INTRODUCTION 

HESS J1858-I-020 is a weak gamma-ray source (1.6% Crab Flux) that was reported not 
to have any clear cataloged counterpart at any wavelength (Aharonian et al. 2008). HESS 
J1858-I-020 is a nearly point-like source, with a slight extension of ~ 5 arcmin along its 
major axis. Its differential spectral index is 2.2 ±0.1. The nearby radio source G35.6-0.4 
was recently re-identified as a SNR (Green 2009). HESS J1858+020 lies towards the southern 
border of this remnant. Paron & Giacani (2010) have found, using the ^^CO (J=l-0) line 
from the Galactic Ring Survey and mid-IR data from the Galactic Legacy Infrared Mid-Plane 
Survey Extraordinaire (GLIMPSE), that there is one or several MCs towards the southern 
border of SNR G35.6-0.4, likely at the same distance of the remnant (10.5 kpc). Paron & 
Giacani (2010) also provide estimates of the cloud's total molecular mass and density. They 
report that the cloud is composed of two molecular clumps, one over the SNR shell and 
the other located at the center of HESS J1858+020 with a molecular mass and a density of 
~ 5 X 10^ M0 and ~ 500 cm~^, respectively. They proposed, using a simplified analytical 
approach described in Torres et al. (2003), that hadronic gamma-ray emission within the 
clouds, produced by protons diffusing away from the SNR G35. 6-0.4, is a possible origin of 
HESS J1858+020. In a more recent paper, Paron et al. (2011) give more details about the 
molecular material, obtained via observations with the Atacama Submillimeter Telescope 
Experiment. They discovered a young stellar object (YSO), probably a high mass protostar, 
embedded in the molecular clump, but no evidence of any molecular outflows which might in 
principle reveal the presence of a thermal jet capable of generating the observed gamma-rays. 
Paron et al. (2011) concluded again that the most probable origin for the TeV gamma-ray 
emission are hadronic interactions between the molecular gas and the cosmic rays accelerated 
by the shock front of SNR G35.6-0.4. Here, we focus on a more in-depth analysis of this 
possibility. 

* dtorres@ieec.uab.es 
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An important ingredient to assess the latter proposal is to have an understanding of 
the GeV emission (if any), of the region. However, the GeV emission at the sky position of 
HESS J1858+020 has not been studied up to now. A look at the Fermi-LAT First Source 
Catalog (Abdo et al. 2010a) shows that there is no source at the position of the SNR or 
the H.E.S.S. source. The Fermi source IFGL J1857.H-0212c is the closest one: it is only at 
0.3° from HESS J1858+020, it is unidentified, has an error radius of 0.08°, and its spectral 
index is 2.31± 0.04. In the IFGL analysis there is no indication that IFGL J1857.1+0212c 
is confused. The pulsar PSR B1855+02 lies close (see the map below), but no gamma-ray 
emission has been detected originating from it in the Fermi-LAT data (Abdo et al. 2010b). 
PSR B1855+02 was suggested as being possibly related with the SNR G35.6-0.4 but the 
association is unclear (see Green 2009 for details). The proximity of HESS J1858+020 and 
IFGL J1857.1+0212c, particularly since the point spread function (PSF) of the Fermi-LAT 
instrument is larger than the separation of the two sky coordinates of interest, requires a 
careful analysis in order to evaluate the possible relations among these sources. 

The case of middle-aged SNRs interacting with MCs, with the subsequent production 
of gamma-ray sources at GeV and/or TeV energies has also been analysed, for instance, for 
W51 (Abdo et al. 2009), W28 (Abdo et al. 2010e, Li & Chen 2010, Ohira et al. 2011), and 
IC 443 (e.g., Abdo et al. 2010f, Torres et al. 2010, see also Rodriguez Marrero et al. 2008). 
However, in these cases there was a GeV detection, except for two of the four TeV hotspots of 
the SNR W28, namely HESS J1800-230A and HESS J1800-230C for which only upper hmits 
were imposed. These instances of GeV - TeV non-simultaneous detections around SNRs are 
probably the most interesting ones from a cosmic-ray viewpoint (see the discussion in Funk 
et al. 2008). They allow us to explore an hadronic origin of the highest-energy gamma-ray 
emission, through cosmic rays accelerated by the SNR shock front, which, diffusing away 
from it, interact with MCs at a certain distance. Given that only the most energetic protons 
reach the separate clouds at a given time, these cases allow us to investigate the diffusion 
environment in which the cosmic-ray propagation proceeds. 

Section 2 is the observational core of the paper and we present there the analysis of more 
than 2 years of Fermi-LAT data of this region. We put special care to analyse whether a 
mis-localization of IFGL J1857.1-|-0212c could render it coincident in projection with HESS 
J 1828-1-020 and prove that this is likely not the case. We then impose upper limits to the 
GeV radiation at the position of HESS J1858-I-020. These limits are used in Section 3 to 
establish possible scenarios where the constraints (detection at TeV, non-detection at GeV) 
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could be satisfied. A discussion of tliese results and a comparison with the similar case of 
SNR W28 with nearby MCs (MCs) is given at the end. 

2 GEV ANALYSIS 

We searched for gamma rays emitted by the HESS J1858+020 analyzing the publicly avail- 
able ~ 28 months of Fermz-LAT data in the energy range 100 MeV - 100 GeV, from 4 
Aug 2008 to 19 November 2010. The search is performed by means of the binned likelihood 
spectral analysis, using the official tool (gtlike) released by the Fermz-LAT collaboration 
[Fermi Science Tools v9rl7p0). All the data, software, and diffuse models used for this anal- 
ysis are available from the Fermi Science Support Center |j Events from the "Pass 6 Diffuse" 
class were selected, i.e. the event class with the greatest purity of gamma rays, having the 
most stringent background rejection (see details in Atwood et al. 2009). The "Pass 6 v3 
Diffuse" instrument response functions (IRFs) were applied in the analysis (Rando 2009). 
We selected events with energy E >100MeV in a square aligned with celestial coordinates, 
inscribed inside a circular region of interest (ROI) of 15° radius, centered on the H.E.S.S. 
source. The good time intervals are defined such that the ROI does not go below the gamma- 
ray-bright Earth limb (defined at 105° from the Zenith angle), and that the source is always 
inside the LAT field of view, namely in a cone angle of 66°. Source detection significance 
is determined using the Test-Statistic (TS) value, TS = — 2(Lo — Li), which compares the 
likelihood ratio of models including an additional source, Li, with the null- hypothesis of 
background only, Lq (Mattox et al. 1996). 

To apply the likelihood analysis, a spectral-spatial model containing diffuse and point- 
like sources was created. Using the IFGL Catalog we have 93 sources closer than 20° and 

5 closer than 3° from the HESS J1858+020 position. For the Galactic diffuse emission we 
used the spectral-spatial model "gll_iem_v02.fit", which is the one used by the Fermi-LAT 
collaboration in order to build the Fermz-LAT First Source Catalog (Abdo et al. 2010a; 
IFGL as referred to hereafter). The isotropic diffuse emission was modelled by the spectrum 
described in the "isotropic_iem_v02.txt" file. The normalization factors of these two models 
were left free in the fitting procedures. The spectral-spatial model also included all the 
point-like sources from the IFGL list closer than 20° to the H.E.S.S. source. Each of those 
point sources was modeled with a simple power-law, with the exceptions of those sources 

^ http://fermi.gsfc.nasa.gov/ssc/ 
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Figure 1. Counts map of Fermi-LAT data for E ^ 1 GeV with subtracted Galactic and isotropic diffuse emissions, in a window 
of 2? 2 centered on HESS J1858+020. A gaussian smoothing with sigma=0.2 is applied. The Fermi-LAT sources from the 
IFGL catalog are marked. The center of the H.E.S.S. sources in the field are labeled with a small cross, while their extension 
is represented by the ellipsoidal fit from Aharonian et al. 2008. The position of PSR B1855+02 is marked with a diamond. 

that are known to be pulsars, for which a power-law with an exponential cut-off was used. 
The spectral parameters of those sources were set at the IFGL values or those from the 
Fermi-LAT First Pulsar Catalog (Abdo et al. 2010b). The flux parameter of all the point- 
like sources closer than 3° to HESS 1858+020 were left free in the likelihood fit, except for 
the closest one, to which we pay special attention below. 

Figure [l] shows the counts map of Fermi-LAT data for E > 1 GeV with subtracted 
Galactic and isotropic diffuse emissions, in a window of 2" x 2" centered on HESS J1858+020. 



2.1 Testing the hypothesis of an association with the closest Fermi source 

HESS J1858+020 is in the Galactic plane at coordinates 1=35.578°, b=-0.58r (RA=284.584°, 
DEC=2.09°, J2000) and does not spatially coincide with any source listed in the IFGL. 
However, this region is crowded with gamma-ray sources. Indeed, in the spectral-spatial 
model created for the likelihood analysis we collected 93 Fermi-LAT sources, with IFGL 
J1857.1+0212C being the closest source to HESS J1858+020. Even though the H.E.S.S. 
source is outside the error box of IFGL J1857.1+0212c, an association can not be excluded 
a priori, because it is known that the Galactic diffuse model used to build the IFGL catalog 
has the highest uncertainties on the Galactic plane, amounting to ^ 6% as estimated in the 
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Table 1. Fermi-hAT upper limits at the position of HESS J1858+020. 



Energy 


^95%UL 






(^95%UL 






^p95%VL 


(GeV) 


10-9x 


ph 


cm" 


-2s-l 


10-"x 


ergs" 


-1 —2 

cm 


lO^-'^^x crgs^^cm^^ 


0.10- 1.00 


50 








1.85 






0.72 


1.00 - 10.00 


6.93 








2.56 






1.00 


10.00 - 100.00 


0.21 








0.76 






0.30 


>0.10 


113 








7.33 






- 



case of the supernova remnants W51C and W49 (Abdo et al. 2009, 2010c; see also §4.7 of 
the IFGL catalog). Furthermore, we are analyzing this region using more than 2 years of 
Fermi-LAT data, against the 11 months used for the IFGL, which may induce changes in 
the position and fluxes of the sources in the spatial-spectral model. 

To investigate the hypothesis of a possible association of IFGL J1857.1+0212c with 
HESS J1858+020, we performed two likelihood analyses. In one we set in the spectral- 
spatial model the coordinates of IFGL J1857.1+0212c at the IFGL position. In the other, 
we set it at the position of the H.E.S.S. source. The analyses are performed setting free the 
flux parameters of all the point-like sources closer than 3° to HESS J 1858+020, and mod- 
eling IFGL J1857.1+0212c with a power-law with index and flux free. The only difference 
between both analyses was the assumed position of IFGL J1857.1+0212c. This method al- 
lows testing whether the position listed in the IFGL Catalog is sustained. In the flrst case we 
obtained a test statistic value TS=863, while in the second case, with IFGL J1857.1+0212c 
displaced, we obtained TS=509. The comparison of the two TS values obtained changing the 
coordinates of IFGL J1857.1+0212c suggests that the association is signiflcantly unlikely, 
given that ATS > 350. 

In order to take into account the systematics due to the uncertainties of the Galactic 
diffuse model, the analyses were repeated by artiflcially changing the normalization of the 
Galactic diffuse model by ±6% (see. e.g., Abdo et al. 2010d). Also in these cases, the dif- 
ferences of the TS values (for the putative change in position of IFGL J1857.1+0212c) was 
large {ATS > 300), and suggests that the association is signiflcantly unlikely. Thus, we be- 
lieve we can safely entertain the hypothesis that HESS J1858+020 and IFGL J1857.1+0212c 
are not associated and proceed to impose upper limits. 

2.2 Upper limits 

Once the hypothesis of association is rejected, the signiflcance of the plausible gamma-ray 
emission from HESS J1858+020 was evaluated by means of adding an extra source at its 
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position in the spectral-spatial model for the likelihood analysis. We model it with a power- 
law. When redoing the analysis, we set free the photon index and flux parameters of HESS 
J1858+020, and the flux parameter of all the point-like sources closer than 3°, with the 
exception of IFGL J1857.1+0212c. As discussed above, this is the closest IFGL source to 
the position of interest; and for it, we flxed all its parameters. The latter choice can be 
understood if one takes into account that the Fermi-LAT PSF at 1 GeV (around the peak 
of the LAT sensitivity) is 0.8°. If a new source is supposed to be at a smaller distance from a 
bright one that is already known in the spectral-spatial model, we expect that the likelihood 
fltting procedure will not suppress it even if fake, but will simply share the photon counts 
among them, resulting in a good signiflcance (TS>25) for both sources. Such case has been 
found, for instance, for SGR 1627—41 (Abdo et al. 2010d), where it was found that the high 
TS derived by the gtlike analysis has been caused by the presence of the rather strong 
unidentifled source (IFGL J1636.4— 4737), which lies at 0.12° from the magnetar (although 
as in this case, it was positionally incompatible with the magnetar). 

With these settings, the HESS J1858-I-020 likelihood analysis always results in a TS <25, 
implying no detection. After fltting, we derived 95% flux upper limits for E >100 MeV with 
the proflle method, increasing the flux obtained for HESS J1858-I-020 until the maximum log 
likelihood decreases by 2.71/2. In the same way, the 95% flux upper limits were evaluated 
in three energy bins (0.1 < E <100 GeV) flxing the photon index to 2.26 (the result of 
the likelihood flt is 2.26 ± 0.13). The uncertainties of the Fermi-LAT effective area, and of 
the Galactic diffuse emission are the two main sources of systematics that can affect the 
results derived with our analysis. We estimated these systematics effects by repeating the 
upper limits analysis using modifled instrument response functions that bracket the "Pass 
6 v3 Diffuse" effective areas, and changing the normalization of the Galactic diffuse model 
artiflcially by ±6%. In addition to the evaluation of 95% flux upper limits (i^95%UL^^ ^^ 
derived the energy flux upper limits (G'95%ul _ j^^E{dN/dE)dE)^ and the SED points 

In the same way, we have also evaluated the upper limits on the base of Helene's Bayesian 
method (Helene 1983). The main difference respect the proflle method is that in this case 
the 95% upper limits are found integrating the likelihood proflle (function of the source flux 
F) starting from F = 0, without any assumption on its distribution. Whereas the results for 
the upper limits obtained with the two methods were similar, the Bayesian one gave those 
more conservative. The results of this analysis are reported in Table 1. 
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2.3 Additional checks 

We additionally tested that by changing the energy threshold (from 100 to 200 MeV), the 
Fermi Science tools (from v9rl7p0 to v9rl8p6), slightly different selection cuts in gtmktime, 
and different number of sources in the background model (including those in a ROI with size 
from 10° to 20°), all results are stable. Additionally, we have repeated the analyses using 
an updated version of the Fermi-LAT Catalog sources, built using two years of data and 
thus more compatible with the data of the HESS J1858+020 region we focus upon. With 
this catalog (which is yet internal to the Fermi-LAT collaboration), the number of sources 
closer than 20° from the HESS J1858+020 source are 107, while those closer than 3° are 
10. No significant change in the previous results where found. We have also considered that 
some of the closest sources could be extended, and also found the results to be stable. We 
have also re-computed the upper limits with different values of the assumed photon index 
and found the results to be stable too. For instance, assuming a value of 2.1 for the photon 
index, which is the mean value of the diffuse emission spectrum, only affects the upper limit 
in the highest energy bin by about 1%. 

3 MODELS 

In this Section we follow the analysis of Aharonian & Atoyan (1996) for a point-like injection, 
and of Li & Chen (2010) for a non-point-like injection of cosmic-rays in order to assess the 
hadronic production of gamma-rays in the proximity of SNR G35.6-0.4. See these references 
for notation and further clarifications. In addition of those cosmic-rays injected by the SNR, 
we also consider the contribution of diffuse Galactic protons, as measured in the Solar 
neighborhood (see, e.g., Dermer 1986). For the gamma-ray spectrum calculations we use the 
analytic photon emissivity dN^/dE^ developed by Kelner et al. (2006). 

In order to apply the non-point like approach we consider that the high-resolution radio 
image of SNR G35.6-0.4 (Green 2010) shows an extent of 15x11 arcmin^ and thus we 
approximate the average radius as -Rsnr ~ 20pc using the distance of 10.5 kpc based on 
the proximity of the remnant with the HII region G35. 5-0.0 (Paron & Giacani 2010). The 
remnant radius and age at the transition from the Sedov phase to the radiative phase are 
given by (Lozinskaya 1991): 

tcooi = 2.7 X 10^E°f no°- V, (1) 
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Figure 2. Point-like injection model for the gamma-ray spectrum of HESS J1858-h020, compared with H.E.S.S. observations 
(Aharonian et al. 2008) and Fermi-h AT upper limits (this work). The parameters we used in this figure are in correspondence 
with those given in Table 2. The top panels show cases with impulsive injection. The bottom panels show cases with continuous 
injection. In both panels, from left to right, we show results with increasing value of x> from 0.001 to 0.1. In each panel, solid, 
dashed, dotted, and dash-dotted curves stand for results with increasing values of separation between the accelerator and the 
cloud, according to Table 2. 

respectively. If we assume the SNR explosion energy -Esnr = lO^^erg and the ambient gas 
density no of order 1 cm~^, G35.6-0.4 seems to be near the transition time, of age ~ 27 kyr. 
Thus we use the Sedov (1959) law for the previous evolution of G35. 6-0.4: 

i?SNR = 0.34f^y'7-y''pc. (2) 

We shall adopt different values for the diffusion coefficient, taking into account the slow 
diffusion correction around SNRs are possible (e.g., Fujita et al. 2009, Torres et al. 2010) 
and an efficiency rj = 0.1 (e.g., Blandford & Eichler 1987), i.e., a total injection in cosmic-rays 
of 10^° erg (10% of the explosion power) along the SNR lifetime. In any case, the transition 
from the Sedov phase to the radiative phase is not sharp; and as usual, key parameters 
(such as distance, accurate radius, ambient density, and explosion energy) are uncertain. 
Because of this, we also explore possible fits with the simple point-like injection approach, 
more appropriate for the larger timescales. 



3.1 Point-like injection 

We first adopt the point-like injection approximation and analyse which are the possible 
combinations of parameters, if any, which would show the GeV-TeV phenomenology seen. 
Figure |2] shows some examples of these cases. In the first place, we show the resulting SED 
in the gamma-ray domain for a point-like impulsive cosmic-ray injection, of 10^*^ erg. The 
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Table 2. Parameters of the point-like diffusion model applied to HESS J1858+020, for an impulsive and a continuous injection 
of 10^" erg along the SNR lifetime. 
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Figure 3. Top Panel; Examples of the cosmic-ray spectrum at the position of HESS J1858-I-020 resulting from the point-like 
models of Figure 1 as compared with the cosmic-ray sea (red horizontal line). From left to right, we plot the cases with 
X = 0.001,0.01 and 0.1 presented in Figure 1 and Table 2. Bottom panel: ihid. but for continuous injection. 
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accelerator is placed at the center of the SNR G35.6-0.4, and the SED is the result of the 
interaction of the diffusing cosmic rays with the molecular material, assumed to be located 
in a cloud at a separate position. Table 2 and the top panels of Figure |2] summarize the 
results. The diffusion coefficient is assumed to be in the form of a power-law in energy 
(see, e.g., Aharonian & Atoyan 1996) D{Ep) = lO^^xl-^p/lOGeV)*^ cm^ s"-*^, where x is the 
correction factor for slow diffusion around the SNR (e.g., Fujita et al. 2009) and 6 ~ 0.3-0.7 
(e.g., Berezinskii 1990). The value lO^'^xcm^s"^ is referred to as Diq below. To simplify, we 
have adopted a fiducial value for the slope of the injected cosmic-ray spectrum p = 2, and 
6 = 0.5, since other values within the corresponding typical phase space of these parameters 
do not change our conclusions. Instead, we show results for different values of Dio, from 
10^^ to 10^^ cm^ s^^. We see that the higher the value of Dio the more important becomes 
the contribution of the cosmic-ray sea with respect to that of the accelerator. This happens 
because in order to have a reasonably good fit to the GeV-TeV data, we need larger values 
of SNR-MC separation and MC masses. The latter cannot attain any needed value, of course, 
but it is constrained to be less than the measured one in the region, which is of the order of 
lO'' M0 (Paron & Giacani 2010). In addition, the separation cannot be arbitrarily large, since, 
unless the clouds are significantly in the foreground or background of the accelerator -which 
admittedly may well be the case- the projected distance would be larger than measured. 
The projected distance of the MC to the SNR center is ~ 30 pc. These constraints exclude 
most models with an impulsive injector and Diq > 10^^ cm^ s~^. A value of Diq = 10^^ 
cm^ s~^ with an impulsive injection from the SNR is untenable, for instance, since it would 
require two orders of magnitude more molecular material than available. For the largest D 
(top panel, right), even considering unrealistic masses and distance, the model fits for large 
separations are above the upper limits and seem inadequate. Reasonably good fits, i.e., 
resulting in separation and MC masses in agreement with measurements, and within the 
assumptions made, would require a very low value of the diffusion coefficient. Such values of 
DiQ, representing a very slow diffusion, have been used for instance by Gabici et al. (2007), 
but it is likely that they require a much denser cloud to be feasible. 

In Figure |2] and Table 2 we also show several examples of the resulting SED for a 
continuous injection of 1.17 xlO'^^ erg s^^ (totalizing 10^° erg in the lifetime of the SNR) 
with the accelerator also located at the SNR center. Again, reasonably good fits can only be 
obtained with low values of Diq, up to ~10^^ cm^ s^^. Values of the needed molecular mass 
in interaction with the cosmic-ray population increases with Dio, but in this case, values up 
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to DiQ = 10^^ cm^ s~^ are able to produce a good fit, whereas higher ones would require 
masses beyond those measured. 

In Figure [3] we show the cosmic-ray spectrum at the position of HESS J1858+020 result- 
ing from the point-like models of Figure 1, as compared with the cosmic-ray sea. From left 
to right, we plot the cases with x = 0.001,0.01 and 0.1, using the same color coding and 
parameters as those presented in Figure 1 and Table 2. The only solutions that may be in 
agreement with the total molecular mass measured in this environment correspond to the 
leftmost top panel, and the two leftmost bottom panels in this figure. The particles injected 
by the accelerator dominate the cosmic-ray sea for energies above ~ 10 GeV in these cases, 
and quickly overcome it by several orders of magnitude for higher energies. 

We note however that the caveat of this analysis might be in its point-like assumption. 
The size of the SNR shell, and position of the cloud with respect to the center of the SNR 
are comparable, and it is thus an approximation to consider this setup. In order to explore 
further this environment, especially for the continuous injection, we apply the cumulative 
diffusion model to explore the range of parameters with which the source HESS J1858+020 
could be associated with the newly identified SNR G35.6-0.4 leaving a non-detection at GeV 
energies when the SNR is no longer assumed to be point-like. 



3.2 Non point-like injection 

The way in which the cosmic-ray spectrum is constructed in the cumulative model is dif- 
ferent from the more direct point-like approach. In the cumulative model, we are adding 
contributions coming from different distances, and thus diffusing differently. This is espe- 
cially important when the MC is indeed close to the SNR shell. In this case the cosmic rays 
entering the cloud arrive from distances that span from a few to tens of pc (for those coming 
from the other side of the shell). 

With these SNR parameters, the GeV-TeV data of HESS J1858+020 are fitted (as 
plotted in Figure 111) and the resulting parameters are listed in Table 3. To compare with 
the former case, we have adopted an initial spectral index of escaping protons, p = 2, 
predicted by the classic Fermi-type acceleration process, and 6 = 0.5. We have also explored 
steeper values oi p to fit the gamma-ray spectrum of this source so as to study the parameter 
dependence. For p- values equal to 2.1 or 2.2, the needed 6 would decrease a little in the range 
0.3-0.5 to better fit the spectral data, while the lower limit of -Re would still be at ~40-60 
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Figure 4. Non point-like injection model for the gamma-ray spectrum of HESS J1858-I-020, compared with H.E.S.S. obser- 
vations (Aharonian et al. 2008) and Fermi-hAT upper limits (this work). From left to right, we show results with increasing 
value of Xi from 0.001 to 0.1. The parameters we used in this figure are in correspondence with those given in Table 3. Solid, 
dashed, dotted, and dash-dotted curves stand for results with increasing values of separation. Bottom panels: 



Table 3. Parameters of the non point-like (3D) diffusion model applied to HESS J1858-I-020, for an impulsive and a continuous 
injection of 10^" erg along the SNR lifetime. 



kind 


Re 

(pc) 


P 


S 


Mel 

(lO^M^) 


X 


3D 


22 


2 


0.5 


1.1 


0.001 


3D 


23 


2 


0.5 


1.1 


0.001 


3D 


24 


2 


0.5 


1.2 


0.001 


3D 


25 


2 


0.5 


1.6 


0.001 


3D 


30 


2 


0.5 


6.7 


0.01 


3D 


35 


2 


0.5 


8.2 


0.01 


3D 


40 


2 


0.5 


12 


0.01 


3D 


45 


2 


0.5 


19 


0.01 


3D 


50 


2 


0.5 


70 


0.1 


3D 


60 


2 


0.5 


80 


0.1 


3D 


70 


2 


0.5 


100 


0.1 


3D 


80 


2 


0.5 


120 


0.1 



pc. Higher values of p, e.g., 2.4 and beyond, would produce S < 0.3, which is unreasonable 
for the cosmic-ray diffusion process. 

The fitting results for x = 0.001, 0.01, and 0.1 are shown in Figure 3. As expected, among 
the fitted parameters, the MC mass (Md) and the distance from the SNR center to the MC 






Figure 5. Examples of the cosmic-ray spectrum at the position of HESS J1858-I-020 resulting from the non point-like models 
of Figure [4] as compared with the cosmic-ray sea. From left to right, we plot the cases with x = 0.001,0.01 and 0.1 presented 
in Figure g] and Table 3. 
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(Rci) decrease when the x value is reduced. In the higher end for the diffusion coefficient 
{x = 0.1), the fitted -Rci value should be in the range of 50 pc to 80 pc. But also for the 3D 
model it can be noted that the MC mass required for x = 0.1 is of order 10^ Mq, which 
significantly conflicts with the observed value. In the more realistic 3D model, also the case 
with Dio = 10^^ cm^ s~^ requires a mass beyond the one measured, by a factor of ~ 6 or 
more if the clouds are more distant. On the other hand, in the lower end of Diq {x = 0.001), 
the distance from the SNR center to the MC is quite similar to the SNR radius itself (~ 20 
pc), implying that the MC is essentially at the end of the shock front, which actually could 
be argued from the projected maps of Paron & Giacani (2010). Thus we find that in order 
for a hadronic origin of the gamma-ray emission seen be viable, the diffusion coefficient near 
this SNR should be greatly suppressed with respect to the average Galactic value, by more 
than one order of magnitude. 

Figure |5] shows the results for the cosmic-ray spectrum. Note especially the differences 
at the lowest value of x? the leftmost panel, when compared with the point-like injection. 
The 3D analyses more realistically describes the obtained data at GeV and TeV energies, 
with a much lower enhancement of cosmic rays in the region; about a factor of 100 beyond 
the cosmic-ray sea level at ~ 1 TeV and a crossing at ~ 10-100 GeV. 



3.2.1 An energy-independent diffusion? 

We checked whether the case of an energy-independent diffusion, i.e. (5 = 0, can provide a 
good fit. To do that we have used the same parameters (except for 5) as those shown in 
Table 3 to explore their matching to the multi-wavelength data. As an example, results for 
X = 0.001 and 0.01 are shown in Figure |6| Additionally, we have also set other parameters 
free (apart from fixing (5 = 0) to fit the H.E.S.S. data and the Fermz-LAT upper-limits: For 
instance, in order to fit the H.E.S.S. data with gamma-ray index ~ 2.2, the proton spectrum 
index p should also be 2.2 if 5 = 0. The right panel of Figure [6] shows the case for i?=30 pc, 
p = 2.17, 5 = 0, Mci = 2 X lO'^ Mq, and x = 0-01, which matches the H.E.S.S. data, at the 
price of unavoidably violating the Fermi-LAT upper-limits. The required MC mass would 
be much larger than the one measured in the region, and the model is already ruled out by 
this. A spectral break between GeV and TeV is needed to fit the spectral data, which can 
be naturally produced by an energy-dependent diffusion model. 
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Figure 6. Examples of the gamma-ray spectrum at the position of HESS J 1858+020 resulting from the non point-like models 
of Table 3 but with an energy-independent diffusion. From left to right, the two first panels show the cases with x = 0.001 and 
0.01. The third panel shows the case for with parameters R=SO pc, p = 2.17, ^ = 0, Md = 2 x 10^ M©, and x = 0.01. 



4 DISCUSSION 

The GeV counterparts of some unidentified very high-energy sources have been searched 
by Fermi-LAT and some of them show spectral breaks between the GeV and the TeV 
band (e.g. HESS J1834-087/HESS J1804-216), which may imply that the spectra cannot be 
treated as a single emission component (Tam et al. 2010). These scenarios have been studied 
earlier by Funk et al. (2008), using EGRET data. The very faint X-ray emission of these 
unidentified sources indicates that, unless evolved pulsar wind nebulae, these sources are 
likely not gamma-ray emitters themselves, and cosmic-ray diffusion origins such as the one 
explored in this paper could be a possible alternative. 

The cumulative diffusion model predicts concave gamma-ray spectra, in which two peaks 
are present. The one at the lowest energies is attributed to the diffuse Galactic protons, while 
the one at the highest energies, and putatively responsible for the emission detected from 
HESS J1858-I-020 is attributed to the accelerated protons that escaped from the nearby SNR 
G35.6-0.4. It can be seen that the high-energy emission peak shifts to higher energies when 
the distance between the SNR and the cloud increases. Therefore, accelerator-illuminated 
gamma-ray sources, which have the property of being TeV bright but GeV faint, can be 
explained by diffusion effects. Such a scenario was applied, for e.g., by Li & Chen (2010) to 
explain the gamma-rays of the four sources in the SNR W28 field. In W28, four gamma-ray 
sources with various GeV and TeV brightnesses are accounted for by assuming different 
distances from the SNR center. HESS J1858-I-020/SNR G35.6-0.4 might in principle be 
another case of such generic phenomenology. 

In this paper we have studied HESS J1858+020 at GeV energies, while considering its 
possible association to the radio source, recently re-identified as a SNR G35.6-0.4. The latter 
is found to be middle-aged (~ 30 kyr) and has nearby MCs (Paron & Giacani 2010). The 
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Fermi-LAT First Source Catalog did not list any source coincident with the position of 
HESS J1858+020, but the closest Ferrm-LAT source, IFGL J1857.1+0212c, could have in 
principle been associated to it. A detailed analysis of 2 years of Fermi-LAT data disfavors 
this association, and we have imposed upper limits to the GeV emission from the region of 
interest. Using both, a point-like and 3D models for the cosmic-ray injection and propagation 
out from the shell of the SNR, we consider whether the interaction between SNR G35.6- 
0.4 and the MCs nearby could give rise to the TeV emission of HESS J1858+020 without 
producing a GeV counterpart. We have found that although the phase space in principle 
allows for such a situation to appear, usual and/or observationally constrained values of 
the parameters (e.g., diffusion coefficients and cloud-SNR shell distance) would disfavor 
it. Specifically, at the currently assumed distance of SNR G35.6-04 and the MC complex, 
the diffusion coefficient near this SNR should be greatly suppressed with respect to the 
average Galactic value, by more than one order of magnitude, in order for the gamma-ray 
phenomenology detected to be a viable outcome of the three hadronic models considered 
here. This may be possible, and such cases have earlier been considered in the literature, 
but the size and density of this particular cloud may be too low to generate such a slow 
diffusion timescale. 
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